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Abstract 

We analyze 2-point functions in the relativistic Bose gas on the lat- 
tice, i.e., a charged scalar <p A field with chemical potential /i. Using a 
generalized worm algorithm we perform a Monte Carlo simulation in a 
dual representation in terms of fluxes where the complex action prob- 
lem is overcome. We explore various aspects of lattice spectroscopy 
at finite density and zero temperature, such as the asymmetry of for- 
ward and backward propagation in time and the transition into the 
condensed phase. It is shown that after a suitable subtraction the ex- 
ponents for forward and backward propagation are independent of \x 
and agree with the mass obtained from the propagator at fj, = 0. This 
holds for [i < /i c and shows that below the condensation transition the 
mass is independent of /i as expected from the Silver Blaze scenario. 
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Introduction 

In recent years several interesting developments for the simulation of lattice 
field theories at finite density were presented (for examples related to this 
work see, e.g., [H [2], [3], SI El El El [8] ) - A key ingredient in these developments 
is a reformulation of the original lattice field theory in terms of new variables 
(dual variables) which usually are loops or surfaces. The dual variables are 
subject to constraints, but can be updated very efficiently with generaliza- 
tions of the Prokof 'ev-Svistunov worm algorithm [S] . Although some of the 
most interesting problems, such as finite density lattice QCD, are not yet 
under control, many interesting new concepts and techniques now enlarge 
our toolbox for such future challenges. 

In this letter we address the problem of spectroscopy in finite density 
lattice field theory. We study 2-point functions for a charged scalar 4 field 
with a chemical potential (relativistic Bose gas). We use a dual representa- 
tion [8] in terms of fluxes where the partition function is a sum over closed 
loops of fluxes with real and positive weights, i.e., the complex action prob- 
lem of the conventional representation is solved in the dual form. Our Monte 
Carlo update for the dual degrees of freedom is a generalization of the worm 
algorithm [9] and was presented in [8]. 

Lattice spectroscopy is based on 2-point functions where suitable mono- 
mials of the lattice fields are separated in Euclidean time. At finite density 
these two-point functions are modified since the chemical potential intro- 
duces an asymmetry in time: Matter which is propagating forward in time 
is enhanced by the chemical potential fi over anti-matter propagating back- 
ward in time. Thus one expects that finite \i introduces an asymmetry in 
2-point functions. Furthermore, in finite density field theory condensation 
transitions can appear when the chemical potential exceeds a critical value 
fj, c (see [U [HJ for an analysis of condensation in the relativistic Bose gas on 
the lattice). In this exploratory study we analyze the propagators in both 
phases, at \i < \i c where we study the asymmetry of the 2-point functions, 
and at \i > \i c where we observe condensation. 



Conventional and dual formulation 

Since later we will discuss also a continuum analysis of the charged scalar 
field, we start with presenting the action in the continuum form, 
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The chemical potential [i couples to the conserved charge of the scalar field, 
but also shifts the mass parameter. This has an important implication for 
the free theory (A = 0): for > m the action is unbounded from below 
and the theory becomes unstable. For finite A, sufficiently large chemical 
potential flips the sign of the quadratic term, giving rise to the Mexican hat 
potential and the transition into the condensed phase. 

In a lattice discretization the continuum action ([I]) assumes the form 

S = EN^| 2 + A I^| 4 -E [e-^titx+v + e M "' 4 ^a-p] ) • (2) 

X \ u=l J 

The first sum is over the sites x of a N% x Nt lattice with periodic boundary 
conditions, the second sum is over the four Euclidean directions v = 1, 2, 3, 4, 
and V denotes the unit vector in ^-direction. The degrees of freedom are 
the complex valued field variables <j> x attached to the sites x of the lattice, 
r/ denotes the combination 8 + m 2 , where m is the bare mass parameter and 
the 4 coupling is denoted by A. All parameters are in units of the lattice 
spacing a, in other words the lattice spacing is set to a = 1 throughout this 
paper. The partition sum Z is obtained by integrating the Boltzmann factor 
e~ s over all field configurations, Z = J D[(p]e~ s . The measure is a product 
over all lattice sites x, with cp x being integrated over the complex plane, i.e., 
D o II,. J ( .r/u,. 2^. 

The lattice partition sum in the conventional form can be mapped exactly 
to a dual representation [8], where the new degrees of freedom are integer 
valued fluxes attached to the links. The dual partition sum is given by 



{k,l} \ x \ v 
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The sum is over all configurations of two sets of flux variables: k x ^ v £ Z 
and l x>u £ INq. The fc-fluxes are subject to a constraint enforced by the 
product over Kronecker deltas (for notational convenience here denoted by 
5 nt o = S(n)). The constraints enforce the conservation of /c-flux at each site 
x, i.e., [kx,v — ^x-v,v\ = f° r an x - The Z-fluxes are unconstrained. The 
configurations of k- and /-fluxes are weighted with various weights, one of 
them including the chemical potential, and the factors W are given by the 
integrals W(n) = drr n+1 e -v r2 -^ r4 ^ w hi c h can be evaluated numerically 
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to arbitrary precision. It is obvious that also at finite [i the weight factors 
are real and positive, and the complex action problem is solved in the dual 
representation. 

For the Monte Carlo update we use a generalized worm algorithm [9j 
that was presented in [8]. In that paper thermodynamical observables such 
as the particle number n, the field expectation value (|</>| 2 ), and the corre- 
sponding susceptibilities were studied. All these observables can be obtained 
as derivatives of In Z with respect to the parameters \x and r\ and in the dual 
representation take the form of expectation values of moments of the flux 
variables. Based on these observables the phase diagram of the relativistic 
Bose gas and condensation phenomena were studied. 



2-point functions in the dual representation 

Using continuum notation the 2-point function needed for spectroscopy is 
given by 

C(t) = [ d 3 y (0(^)0(0, 0)*> « e~ Et , (4) 



where we project to zero spatial momentum, such that the energy E of the 
particle at rest, i.e., its mass, determines the exponential decay of C{t). 
Thus, on the lattice we need to evaluate 2-point functions which in the 
conventional lattice representation are given by 



If 1 

'//"-./ - y / D[<f>]e~ ip y <f>* z = - , . (5) 



Here Z y>z denotes the lattice partition sum with two field insertions 4> y (fil- 
It is straightforward to generalize the derivation in [8] to find the dual rep- 



resentation of Z VjZ , 



J y,z — 
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(6) 



At the positions y and z where the source fields <\> y (j>* are located the index 
of the weight factors W(n) is increased (see the Kronecker deltas in the 
argument of the W). Also the constraints for the k- fluxes change to 

II *M E \^ x ' v ~ kx-v,v\ - $x,y + $x,z ■ (7) 
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In the presence of the source fields <f> y 4>* z the usual local flux conservation 
Ylv [kx,v— ^x-u,v\ = at all sites x is modified such that the site z constitutes 
a source of flux, the site y a sink. Thus in the dual representation for Z y z the 
set of configurations of allowed /c-flux consists of closed loops and a single 
open line of flux connecting the sites y and z. 

To evaluate Z y z efficiently we follow the strategy of [5] and define a 
generalized partition function Z where we sum over all possible positions of 
the field insertions, 

Z = ^u,v ■ (8) 

u,v 

The configurations that constitute Z consist of closed loops of flux plus a 
single flux line with open ends at arbitrary sites u, v of the lattice. For the 
Monte Carlo simulation of Z we can thus re-use the worm algorithm of [8], 
where now every step of the worm constitutes an admissible configuration 
of Z (not only the configurations where the worm has closed that contribute 
to Z). The two point functions are then obtained as 

K<PM ~ Z (S u , v W(f u )/W(f u + 2)) z ' W 

where (..)z denotes the expectation value with respect to the enlarged en- 
semble Z. By f u = Y. u [\ k u,v\+\ k u-v,v\ ~b 2(Z Ujl/ + lu—v,v)\ we denote the com- 
bined k- and /-flux at site u that enters the weights W, and the reweighting 
with W(f u )/W(f u + 2) in Q is necessary to correctly obtain the original 
partition sum Z in the denominator of Q. When taking into account the 
projection to zero momentum we end up with the following dual represen- 
tation of the correlator C{t): 

fl(f\ _ ( $t,U4-V4 )Z ,-.„n 

1 J " (5u,vW(f u )/W(f u + 2)) z ■ { ) 

In our simulations we use x Nt lattices with N s ranging from 16 to 32, 
and Nt from 32 to 100. A full update sweep is the combination of one local 
Metropolis sweep for the l xu variables and one worm for the k XfV fluxes. We 
typically use between 5 x 10 5 and 10 6 measurements (10 6 in the plots we 
show) separated by 5 full update sweeps. The number of equilibration sweeps 
is typically 2.5 x 10 4 . All errors we quote are statistical errors determined 
with the jackknife method. The weig hts W in ^ were evaluated numerically 
using Mathematica and pre-stored for the Monte Carlo simulation. 
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Analysis of the free case 



We begin with an analysis of the free case, i.e., we set A = 0. In this case 
we can compute the correlators C(t) exactly using Fourier transformation. 
The purpose of studying the free case is twofold: The exact analytic result 
allows us to test the dual approach and the correct implementation of the 
two point function, but also serves to understand the general behavior of 
two point functions at non zero chemical potential. That latter aspect is 
most transparent in the continuum which we discuss explicitly (for the com- 
parison to the results from the dual simulation we will of course use Fourier 
transformation on the lattice). 

A straightforward calculation gives the continuum correlator C(t) as the 
following Fourier integral over the 4-component of the momentum, 



C{t) 



dp 4 



a ip4t 



(11) 



Factorizing the denominator one can make explicit the poles in the complex 
P4-plane: p\+m 2 — fi 2 + i2fip^ = [p4 — i{m— /i)][p4+i(m+/z)]. Thus we have 
two poles on the imaginary axis at i(m — fx) and at —i{m + fi). For vanishing 
chemical potential the poles are equidistant from the origin at ±im. Finite 
chemical potential fi shifts both poles downwards on the imaginary axis. 



The integral (11) can easily be computed using the residue theorem. For 



forward propagation, i.e., for t > one closes the integration contour in the 
upper half of the complex plane and finds C(t) = exp(— (m — fi)t)/2m, while 
for negative time t < one finds C(t) = exp((m + fj,)t)/2m. For finite T the 
temporal extent becomes finite with length /3 = 1/T. The corresponding 
propagator is obtained by making the f3 = oo propagator periodic, 



1 

2m 



-(m—jxjt 



1 



-(m-/i)/3 



+ 



_ e -(m+n)P 



, (12) 



where the last expression was obtained by inserting the above results for for- 
ward and backward propagation into the sum and summing up the emerging 
geometric series. Thus for finite temperature the propagator is a superposi- 
tion of a forward propagating part with an exponent of m — fj, and a piece 
running backward in time with exponent m + (jl. 



Note that the result (12) is valid only for fj, < m, since for \x > m the 



action ([T]) at A = is unbounded from below and the path integral does not 
exist. For the poles in the complex plane this instability corresponds to a 
transition of the pole i(m — fi) into the lower half of the complex plane. 
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Figure 1: The normalized propagator C(t)/C(0) for the free case at rj = 9.0 
(i.e., m = 1.0) on 32 3 x 100 lattices. We show the results at different 
values of \x and compare the Monte Carlo data (circles) to the propagator 
as computed from (lattice) Fourier transformation (crosses). 

Let us now come to the discussion of our numerical results for the free 
propagator and their comparison to lattice Fourier transformation. We are 
interested in studying the situation at low temperature and use Nt = 100 
where the temperature in lattice units is T = 1/Nt = 0.01, much smaller 
than any other scale in the problem. The spatial volume is 32 3 and we work 
at A = 0, 7] = 9, i.e., the mass is m = 1.0. This is the upper limit of the 
chemical potential \x we can consider before the theory becomes unstable. In 
Fig. [T] we show the propagator C(t)/C(0) (i.e., normalized to 1 at t = 0) as a 
function of t for different values of fj,. Circles are used for the results from the 
dual simulation, and crosses for the data from lattice Fourier transformation. 

The propagators show the asymmetry between forward and backward 
propagation as expected from the continuum case we discussed. It is obvious 
that the dual Monte Carlo data and the Fourier results agree very well. This 
shows that the implementation of the propagator evaluation based on the 
generalized worm algorithm is correct and that we can reliably determine 
2-point functions also at finite chemical potential. 
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Figure 2: The normalized propagator C(i)/C(0) at rj = 7.44 and A = 1.0 on 
lattices of size 32 3 x 100 for different values of \i. 



Non zero quartic coupling 

We now turn on the quartic coupling. In this case the path integral exists 
also for large values of the chemical potential and we can study the change 
of the propagators across the transition into the condensed phase. We work 
on lattices of size 32 3 x 100 and present results for A = 1.0 and rj = 7.44. 
For these parameters the condensation transition was found [8] at a critical 
chemical potential of jx c = 0.170(1). For ji < ji c thermodynamical observ- 
ables were found |8] to be independent of fj,, a property sometimes referred 
to as Silver Blaze behavior |10| . 

In Fig. [2] we show the normalized propagator C(i)/C(0) as a function of 
t and compare the results for different values of the chemical potential ji. 
Again, as expected, the propagator is symmetrical for ji = 0, while for ji > 
we observe a difference between forward and backward propagation. As in 
the free case we find very clean exponential decay showing up as straight 
lines due to the logarithmic scale used on the vertical axis. For forward 
propagation the slope decreases with increasing fj,, while it increases for 
the backward running part. The slope for forward propagation vanishes as 
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0.1685(2) 
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0.15 


0.0196(3) 




0.1696(3) 





Table 1: Fit results for the exponents E + and E- for forward and backward 
propagation at different values of the chemical potential fi at r] = 7.44 and 
A = 1.0 on lattices of size 32 3 x 100. In the last two columns we add 
(subtract) the chemical potential \i from the fit values E + and E-. 



the chemical potential reaches the critical value [i c where the propagator 
becomes constant. For values of fj, above fi c , i.e., in the condensed phase, 
the propagator remains constant. 

Next we analyze if the slopes of forward and backward propagation can 
be understood quantitatively. We determine the exponents E + and EL for 
forward and backward propagation from fits according to Q for several 
values of the chemical potential /i. The results are listed in the second and 
third column of Table [T] According to the Silver Blaze scenario |10| one 
expects that the mass of the excitation that C{t) couples to is invariant for 
different values of /i < fi c . To check this expectation we add fi to E + and 
subtract it from E-. The results for these reduced exponents at different 
fi < n c should be constant in fi and agree with the mass obtained from 
the propagator at fi = 0. The fourth and last columns in Table [T] show 
that this is indeed the case: The numbers for the reduced exponents are all 
close to 0.169, thus coinciding with the exponent of the propagator at \i = 
and agree well with the critical value fj, c = 0.170(1) determined in [8j. We 
conclude that the propagators are understood also quantitatively and the 
mass of the lowest excitation indeed shows the Silver Blaze property, i.e., it 
is independent of \x up to fi c . 

Let us finally comment in more detail on the behavior in the condensed 
phase: There the correlator should be dominated by the Goldstone mode 
which gives rise to algebraic decay, while the massive mode decays exponen- 
tially. However, at finite volume there is no spontaneous symmetry breaking 
and a finite volume Monte Carlo simulation only produces an "erratic con- 
tinuous drift through the degenerate vacuum states" [TTj. Two successful 
solutions to the problem are to either introduce a small external field to 
explicitly break the symmetry [12] or to use interpolators in the 2-point 
functions where <f> x is projected onto the dominant direction of the magne- 
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tization M = Y2 X f° r eacn configuration individually [TT]. In both cases 
one needs the magnetization M. The problem with the dual representation 
we currently use is that we can only access observables that are even in the 
field (f>, while M is odd. A more general dual representation of the model 
with magnetic field is possible where additional sets of dual variables, so 
called monomers, are introduced. The dual configurations are open strings 
with monomers at their endpoints. We hope that with this generalized 
dual representation we can successfully address the problem of studying the 
massless and massive modes of the condensed phase in the near future. 

Summary and discussion 

In this work we have analyzed 2-point functions for a charged scalar 4 field 
at finite density using the lattice formulation. The complex action problem 
of the standard form was overcome by mapping the partition sum to a 
flux representation where only real and non-negative terms appear. The 
2-point functions were mapped to the dual formulation and we evaluated 
them at low temperatures and finite \x in a Monte Carlo simulation based 
on a generalized worm algorithm. 

In a first step we studied the free case to check the correctness and accu- 
racy of the dual approach to 2-point functions. We observe the asymmetry 
between forward and backward propagation as expected from the analysis 
of the free case in the continuum and find perfect agreement of the dual 
simulation with the result from lattice Fourier transformation. 

For the case of non zero quartic coupling we again find the asymmetry of 
the propagators. As n approaches its critical value \x the exponential decay 
of forward propagation weakens and vanishes at the condensation transition. 
In the condensed phase the propagator becomes constant. We find that after 
adding (subtracting) \i from the exponent E + (E-) for forward (backward) 
propagation we get a number independent of [i which coincides with the 
mass obtained from the propagator at \i = 0. This holds for all \i < \i c and 
shows that in this region the mass is independent of \x as expected from the 
Silver Blaze scenario. 

The analysis presented in this study of course only has an exploratory 
character, but we expect that the findings (and maybe some of the tech- 
niques) are relevant also for analyzing other systems with a Silver Blaze 
behavior, such as QCD at finite \x and low temperatures. There the chem- 
ical potential couples to quarks and the lowest excitation is the nucleon. 
Repeating the analysis of the free case for fermions, one again finds the dif- 
ference between forward and backward propagation with exponents m db /i, 
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and it is expected that also in the coupled case after a suitable subtraction 
a mass independent of \x is obtained. While an analysis of this behavior 
is out of reach in full finite density lattice QCD, it would be interesting to 
check for this behavior of nucleon propagators in the QCD related model [3] 
obtained from hopping and strong coupling expansion. 
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